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ABSTRACT 

We calculate the angular momentum transport by gravito-inertial-Alfven 
j^ waves and show that, so long as prograde and retrograde gravity waves are 

O ' excited to roughly the same amplitude, the sign of angular momentum deposit 

in the radiative interior of the Sun is such as to lead to an exponential growth of 
^ ■ any existing small radial gradient of rotation velocity just below the convection 

zone. This leads to formation of a strong thin shear layer (of thickness about 
>- ■ 0.3% R0) near the top of the radiative zone of the Sun on a time-scale of order 

20 years. When the magnitude of differential rotation across this layer reaches 



0^ 



O 



en ■ about 0.1 /iHz, the layer becomes unstable to shear instability and undergoes 

Q ■ mixing, and the excess angular momentum deposited in the layer is returned to 

^ ■ the convection zone. The strong shear in this layer generates toroidal magnetic 

■^ ' field which is also deposited in the convection zone when the layer becomes 

Q_^• unstable. This could possibly start a new magnetic activity cycle seen at the 

O ■ surface. 
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1. Introduction 



The Sun is currently losing angular momentum at its surface via a wind at a rate of 
the order of 10^^ g cm^ s^^, which is slowing down the surface layers. However, we know 
from the frequency splitting of p-modes that the radiative interior of the Sun is rotating as 
a solid body at a rate that is not very different from the mean surface rotation rate. Thus 
the loss of angular momentum at the surface is communicated to the rest of the Sun on a 
time-scale of the order of the solar age or less. 

The convective motions in the outer third of the Sun can very efficiently redistribute 
angular momentum on a short convective time-scale of about a month. However, processes 
responsible for the redistribution of angular momentum in the radiative interior are less 
certain. One possibility is that magnetic torques can extract angular momentum from the 
radiative interior. Charbonneau & MacGregor (1993) have carried out a detailed numerical 
analysis of this transport process. However, there are some drawbacks to that mechanism, 
which are discussed in Zahn (1997). 

Another possibility is that the gravity waves generated in the convection zone can 
extract /deposit angular momentum in the radiative interior. The angular momentum flux 
carried by these waves is enormous and they are very efficient in redistributing angular 
momentum in the radiative interior (Schatzman 1993a-b; Kumar & Quataert 1997; Zahn, 
Talon & Matias 1997). 

Here, we explore this second mechanism in some detail. In particular, we discuss the 
effects of the Coriolis force and of magnetic fields on the dispersion relation of gravity waves 
and on the deposition of angular momentum in the radiative interior on the Sun (§2). We 
describe how these waves lead to the formation of a strong and thin shear layer just below 
the convection zone which becomes unstable when the gradient exceeds a critical value 
leading to mixing of elements and angular momentum with the convection zone (§3). We 
further discuss how this shear layer could contribute to the generation of magnetic field and 
perhaps to the magnetic cycle observed in the Sun (§4). 



2. Angular momentum transport by waves 

We derive below a general expression for the angular momentum fiux in waves. The 
dispersion relation for gravito-inertial waves in the presence of a magnetic field is calculated 
in 52.2. 
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2.1. Angular momentum flux in waves 

Let us consider a fluid Lagrangian density L wliich is a function of the displacement 
fleld I, and of its temporal and spatial derivatives. If L is not an explicit function of the 
azimuthal coordinate 0, the z-component of the angular momentum is conserved (this is a 
particular case of the Noether's theorem cf. Quigg, 1983). The Lagrangian density thus 
obeys the relation 
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It follows from the Euler-Lagrange equation that the last term is zero leading to 
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This equation expresses the conservation of wave angular momentum. The angular 
momentum density and flux are given by: 
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We show below that the above expression for angular momentum density {Cam) has the 
correct sign and magnitude. 

Consider a plane wave in a homogeneous medium traveling along +x direction, 

I = C,qX sin(ti;t — kx). (6) 

The only term in the Lagrangian containing d^/dt is the kinetic energy term, so it is 
straightforward to calculate the angular momentum density for the wave using eq. (H): 



£. 
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at 00 
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dt, 



(7) 



where y = r sin (p. Since the wave is traveling along the +x direction, its linear momentum 
density is also along the +x direction and has the magnitude p = p{uj^o)^{k/u!). The 
angular momentum density of the wave should be Cam = x x p = —pu'^{k/uj)yz, which is 
identical to the expression (0). 
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The wave displacement field in the spherical Sun can be written as 

/ d 3 d \ 

I = Kr, ^hT^, —^TTT yemid, 0) COs{ujt - TrKJ)). 

\ oO sm ad) 



Substituting this in the expression for Cam (eq- H), we obtain 



Integrating Cam over 9 and we obtain 



2c'2 



sin2(cut - m^) + ^^^V^^ cos^(cjt - m(\)) \ . (9) 



(C^) = ^[C^ + ^(^ + i)e^]=— , (10) 

where E is the energy density in the wave. Prograde waves [m > 0) thus carry positive 
angular momentum whereas retrograde waves {m < 0) carry negative angular momentum^. 

Waves deposit their angular momentum in the star only when and where they are 
damped. For gravity waves, the main damping process is the photon diffusion, which is 
discussed in §2.4. The angular momentum flux crossing a spherical shell may then be 
expressed as a function of depth: 

Cam{^J,r) = Cami^J,rc) exp[-r(u;, £, r)] , (11) 

where T{ij,i,r) is the damping "optical depth" for a wave of frequency u, and degree i, 
between the radius r and Vc (the radius of the bottom of the convection zone) cf. §2.4. 



2.2. Dispersion relation for Gravito-Inertial-Alfven waves 

The linearized momentum equation that includes both Coriolis and Lorentz forces is 

P^ + ^P^ - P^S + '2pnxv = ^^1^ X B (12) 

where ft is the rotational frequency, g = —gr is the gravitational acceleration, ^ and v are 
the fluid displacement and velocity associated with the wave, and pi, pi and Bi are the 
Eulerian perturbations of density, pressure and magnetic field associated with the wave. 
We assume for simplicity that the unperturbed magnetic field B and the rotational speed 



^The sign of the angular momentum carried by the gravity waves was improperly defined in Kumar & 
Quataert (f997) and in Zahn, Talon & Matias (1997), as was pointed out by Ringot (1998) 



-5- 



O are constant. In the short wavelength hmit, we can take the spatial and the temporal 
dependence of all of the perturbed quantities to be proportional to exp(zk ■ x — iut). 
Substituting this in the above equation we obtain 



pco'^C, + ik.pi — Pig — 2iujp n, X ^ 



i{kx Bi] 
47r 



X B. 



Making use of the linearized equation of entropy conservation 

Pi 



§ + ^iV^, 



c" g 



(13) 



(14) 



where C,r is the radial component of the displacement vector, and taking k ^ g/c^, i.e. the 
wavelength is much smaller than the density scale height, we obtain 



o , o ^ ^ (k X Bi) _ 
- puj^^ + ikpi + pNXr'c - 2iujpd X I = ^ x B. 



(15) 



Substituting for Bi from the linearized induction equation for a perfectly conducting fluid, 
and making use of the incompressibility condition, k ■ | = 0, we obtain 



- pu'^^ + ikpi + pN^^ri- - 2iujp (fi x |) 



(k-B) 
47r 



{(k-B)|-(^B)k} = 0. 



The vector cross product of the above equation with k yields 

(u^ - (k ■ B) Vj) (k X I) - N^^r (k X r) - 2tuj (k ■ rj) ^ = 0, 



(16) 



(17) 



where V| = B'^/Airp is the Alfven velocity, and B = B/B. The vector cross product of eq. 
(P^) with k gives k x ^ in terms of |, which, when substituted back into eq. (0), results in 



(o;2-(k-B)Vi)'A;2 
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The r-component of the above equation yields the desired dispersion relation: 
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where k = k/k. This equation can be solved easily to determine the wave frequency 
1 



uj' = (k.B)Vi+^ 



2 . . 1 1/2' 



(Nxk)'+4(k-n)'± (Nxk)' + 4(k-0)' +16(k-B)V^(k-ri 



(20) 
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This dispersion relation has two branches. The lower frequency branch corresponds 
to Alfven waves and the high frequency branch, to gravito-inertial-Alfven waves. The 
dispersion relation for the high frequency branch can be written as 

cu' ^ (k ■ B) Vj + (N X k)2 + 4(k ■ 0)2. (21) 

The radial wave number kr, when the magnetic field lies in the horizontal plane, can 
be obtained from the dispersion relation (eq. |T9|) and is given by 



2uj^n^ke sin 26 + [Auj^Vt^k'^g sin^ 26 + (cjf - Auj^Q^ cos^ 6) {{N^ - ujDujfkl + Au'^n'^kl sin^ 6] 



1/2 



where 



Ujf - 4cu2fi2 cos2 6 

(22) 

^1 = ^"^ -{k-BfVl, (23) 

Q, = r2(cos 6r — sin 6 6), kh is the horizontal wavenumber, kg = k^ — fc^, k^ = m/r sin 6, and 
m is the azimuthal number. Figure 1 shows k^ as a function of uj for a few cases of 6, B and 
kh- 

The frequency uj in the expression for the wavenumber is the frequency as seen in the 
local rest frame of the fluid i.e. uo = uji — rriQ, where uj is the wave frequency in an inertial 
frame. The component of the wave- vector along the (p direction being conserved along the 
ray path, the prograde waves moving inwards in a region of increasing rotation speed are 
Doppler shifted to lower frequencies. From the above equation, we see that the wavenumber 
of these waves increases with decreasing frequency (see also fig. 1) therefore enhancing their 
damping. Thus, in a region of increasing rotation speed, prograde waves are dissipated 
more strongly than retrograde waves. 

For pure gravity waves, Zahn et al. (1997) discussed the existence of a critical layer 
determined by the condition uj = m6Q, where 6Q is the difference between the rotation rate 
at radius r and the base of the convection zone. In that case, a wave is completely absorbed 
since its radial wave number diverges and its group speed goes to zero. The presence 
of non-zero magnetic field and/or rotation modifies the condition for the appearance of 
critical layers, which can occur even when tu > in the local rest frame. However, it can be 
shown that uf = ±2ujQcos6, for which the denominator in equation (^) is zero, does not 
correspond to a critical layer since the numerator at this frequency is zero as well, unless 
either Q cos 6* = or the wave frequency in the local rest frame of the fluid {uj — mSQ) is 
zero. For the low frequency waves (which carry most of the angular momentum) radiative 



^Barnes et al. (1998) discuss critical layers for gravito-Alfven waves and their results can be obtained 
from ours by setting = 0. 
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damping is very strong and the waves are damped on a distance scale much smaller than 
the solar radius even when waves do not encounter a critical layer. This is discussed in 
more detail in §2.4. 

Wave propagation requires uj > (k- B) Va- For horizontal magnetic field of strength less 
than about (2-10^/£) Gauss wave propagation in the radial direction requires the frequency 
to be greater than 0.5 /xHz, while for a radial magnetic field of strength less than about 10^ 
Gauss the condition for wave propagation is a; > O.S/iHz (angular momentum fiux in waves 
peaks at the convective frequency of 0.3 /iHz and falls off as z/~^-^ at higher frequencies). 
Waves will undergo reflection when the radial component of the magnetic fleld is sufficiently 
large. 



2.3. Excitation of waves by turbulent convection 

Gravity waves with periods of about 10 days are expected to be excited in the Sun by 
a number of different processes such as the Reynold's stress in the convection zone, plumes 
penetrating in the stably stratified radiative interior, etc. In the Earth's atmosphere, the 
latter process is known to be the most efficient and it might be the dominant process in the 
Sun as well. However, the energy fiux in gravity waves resulting from this process is subject 
to great uncertainty because we know little about the properties of plumes at the base of 
the solar convection zone. We will therefore consider wave excitation only by Reynold's 
stress. The results of this paper (the formation of shear layer at the top of the radiative 
interior of the Sun by gravity waves) can be easily modified to incorporate more efficient 
gravity waves generation by plumes or some other process. Generally, larger fiux in gravity 
waves means more efficient formation of shear layer on a smaller time-scale. 

The energy fiux per unit frequency in waves just below the convective envelope of the 
Sun, as a result of excitation by the Reynold's stress, can be calculated using the following 
expression (cf. Goldreich et al. 1994): 



^^ ~4nJ r^ 




exp 



-""^(^ +''/*'] iri^' <'*> 



where ^r and [i{i + l)]^/^^/i are the radial and horizontal displacement wavefunctions which 
are normalized to unit energy fiux just below the convection zone, v is the convective 
velocity, L is the radial size of an energy bearing turbulent eddy, tl ~ L/v is the 
characteristic convective time, and h^ is the radial size of the largest eddy at r with 
characteristic frequency of u or greater {h^ = Lmin{l, {2ujtl)~^'^}). The gravity waves are 
evanescent in the convection zone, the region where they are excited. The above equation 



was derived under the assumption that the turbulence spectrum is Kolmogorov and, as 
mentioned earher, it ignores wave excitation resulting from convective overshooting. 

The energy flux spectrum just below the convection zone is calculated numerically 
using equation (p^ ) and the result is shown in flg. 2. Also shown in flg. 2 is the frequency 
integrated wave action as a function of wave degree £; angular momentum flux for waves of 
a given m is equal to m times the wave action. We provide below a rough estimate of the 
energy flux in gravity waves in order to gain some understanding of the numerical results. 

For the low frequency waves of interest here, the WKB solution to the wave equation is 
quite accurate and can be used to calculate the energy flux. The displacement wavefunctions 
^r and ^h in the short wavelength limit are proportional to k~^^'^ exp{i J drkr), and the 
energy flux in waves is puj'^{C,r + ^(^ + ^)^h)'^gy "where kr is radial wave number given by 
equation (]22|), and Vg is the wave group velocity. 

Using the continuity of the radial velocity at the interface of the radiation and the 
convection zone we flnd that the properly normalized radial displacement wave function in 

the convection zone is 

ikHN.V/' 



^-^ .;(piViV,)V^ ^"Pi-i^^H> (25) 

where A''^ ~ (g/HY^"^ is the Brunt- Vaisala frequency at the top of the radiation zone and 
Nc is the convective frequency at the bottom of the convection zone. Substituting this 
into equation ( p^ and making use of L ~ if and pv^ ~ Fq (the energy flux carried by 
convection), we flnd the energy flux in waves to be 

fP(.) ~ fe^ig>^P(-g^^i). (26) 

r^Nt 

The subscript u on H emphasizes that it is the scale height at a location in the convection 
zone where the characteristic time-scale of energy bearing eddies is uj~^. For a polytropic 
atmosphere of index n, the density scales with depth z (measured from the surface) as z^, 
and therefore the convective velocity scales as f oc 2;^"/^. Thus the convective frequency 
uj oc 2;^("+^)/^, and the scale height H^ ~ Hc{Nc/oj)^^^"'^^''; He is the scale height at the 
bottom of the convection zone. Substituting this into the above equation we obtain 

,(,), ., F^klHl eM-Hlkl) ( u ^-i5/(n+3) 



Fr(^)^ — — — (^^) . (27) 



" c 



Since the polytropic index of the convection zone is approximately 1.5, we flnd the frequency 
dependence of Fj^ to be ca;"^'^. The numerical calculation of the energy energy flux, using 
equation (|2^) , gives the frequency dependence to be tu"'^'^. The difference arises because of 
our neglect of the frequency dependence of k^ in equation (^) for ^r- 
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(c) 

Integrating F^ over the wave frequency and the horizontal wave number, we find the 
total energy fiux in waves to be given byQ 

4^) -^- FeMu (28) 

where A^f ~ 3 x 10~^ is the turbulent Mach number at the bottom of the convection zone. 
Thus the total energy fiux in the long period gravity waves just below the solar convection 
zone is about 0.03% of the solar luminosity; more accurate numerical calculation based on 
the evaluation of eq. (24) gives the energy fiux to be 0.1% of the solar luminosity. 



2.4. Radiative damping of waves 

For the long period or short wavelength waves which dominate the transport of 
angular momentum, the radiative diffusion of photons from regions of positive and negative 
temperature fluctuation across a wavelength is the dominant damping process. The 
damping opacity between radii r and r^ (the base of the convection zone) is calculated 
using: 

where 

ui^\r)=u + m[Q,-Q{r)], (30) 

is the wave frequency in the local rest frame of the fluid at radius r, Q^ is the angular 
rotation speed at the base of the convection zone, the damping rate 

j{uj,i,m,r)^ Kk^ (31) 

A- = iM! « ?£.^, (32) 

Sp^KCp 5p 

is the thermal diffusivity, a is the Stefan-Boltzmann constant, k, is opacity per unit mass, 
Cp is the speciflc heat per unit mass, F^. is the radiative energy flux, Ht is the temperature 
scale height, p is the gas pressure, and kr is the radial wavenumber of the wave. 

For gravity waves, kr ~ khN/u and Vg ~ uj"^ /{k^N). Therefore, their damping increases 
rapidly with decreasing frequency and increasing £ as uj~'^C.'^. The wave damping length 
dw ~ "^9/7 ~ 5puj/{2FrHTkr), and so the ratio of the damping length and the wavelength is 



T{u;J,m,r)=l dr' ' \ \' [^ \ [^ (29) 



Puj^Ht / uj \^ / r, ^2 
Frkl ^ [Wj [h^ 



^-^^ ^ -^ ~ [wJ ^ ^ ^ ^ ' ^''^ 



•^Goldreich and Kumar (1990) had obtained the same resuh for the energy flux in surface gravity waves. 
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where Nc is the convective frequency at the base of the convection zone. We have used the 
relations Fr ~ p{NcHtY ^-nd A^^ ~ g/ Ht in deriving the above scahng. At the bottom 
of the solar convection zone rd Ht ~ 10, therefore for dipole waves of frequency N^. the 
damping length is about 15 wavelengths. For waves of frequency 0.6 yuHz {uo/Nf. k, 4) the 
wave damping length is of order the wavelength for £ ~ 20 (see fig. 3 for a more accurate 
numerical result). 

Numerical calculation of damping length as a function of wave frequency and d, are 
shown in figure 3. Note that, even in the absence of a critical layer wave damping is 
extremely efficient (see also Press, 1981). Waves with frequency of 0.5/iHz and i = 15 
are damped on a length scale of about 10 wavelengths. Our numerical calculations for 
the evolution of shear layer below the convection zone, presented in §3, include radiative 
damping of waves throughout the shear layer and not just at the critical surface. 

Raditive damping only affects the gravity wave part of the gravito-Alfven-inertial 
waves. Thus, after traversing some distance, a gravito-Alfven-inertial wave can be converted 
into an Alfven-inertial wave as a result of radiative damping. We have however ignored 
this complication here because, for the waves of frequency greater than about 0.6 /iHz 
considered in this paper, the buoyancy force is more important than magnetic tension and 
the Coriolis forces as long as the toroidal magnetic field below the convection zone is less 
than 10^ Gauss, and the latitute is less than 45° (see §2.2). Thus, much of the angular 
momentum fiux carried by gravito-Alfven-inertial waves gets deposited in the fiuid before 
they turn into Alfven-inertial waves. Moreover, it can be shown that for parameters of 
interest to the solar problem, the finite plasma conductivity damps out Alfven-inertial 
waves near the critical surface. 



3. The formation and evolution of a shear layer 

Long period gravity waves generated at the base of the solar convection zone are 
efficient at redistributing the angular momentum in the radiative interior (Zahn, Talon & 
Matias, 1997, Kumar & Quataert 1997). The angular momentum deposited by gravity 
waves just below the convection zone is such as to enhance any pre-existing gradient of 
angular velocity. The physical reason for this, although not intuitive, is easy to understand. 
Consider the case where the angular speed increases with depth below the convection zone. 
The frequencies of prograde waves moving inward in the radiative interior are decreased 
as seen in the local rest frame of the rotating fiuid whereas the frequencies of retrograde 
waves are increased. Thus the damping of prograde waves is enhanced relative to retrograde 
waves and, as a result, positive angular momentum is deposited just below the convection 
zone enhancing the magnitude of the existing differential rotation. This enhancement of 
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angular momentuiii gradient is confined to a thin layer whose thickness is of the order of the 
wave damping length. Initially, the angular speed increases exponentially with time, and 
when the gradient becomes sufficiently large so that much of the angular momentum fiux 
of prograde waves is absorbed in the shear layer below the convection zone, further growth 
becomes linear in time. Retrograde waves, which were only partially absorbed in this layer, 
continue to propagate deeper in the radiative interior and deposit their negative angular 
momentum in a thin layer which is spun down and develops a strong differential rotation. 

The evolution of angular momentum and formation of the shear layer are determined 
by solving the following equation: 



pr 



,dQ dFr 



dt dr 
where F^ is the angular momentum fiux associated with waves 



(34) 



E 



du;^^ 



W, 



u;,i) 



U) 



exp [— r(c<j,£, m)] . 



(35) 



To estimate the time-scale for the initial exponential growth of the angular speed, we 
expand the damping optical depth r 



T{u,i,m) 



dr'y{r,u—m6fl, 



Substituting this back into equation (|35|) we obtain 



d'y 

dr'y{r,u},£)— I drmSVL—— = ro(r, cc;,£)-|-mri(r, cj, 

duj 

(36) 



EJduJ 



uo dri 



smh{i + l/2)ri 



sinh(ri/2) 



exp (-To) = Y,jd^- 



^(c) 



nFlii + i/2l 

3uj 



exp(-ro). 



(37) 

The last equality was obtained by assuming that iri -C 1 and £ :^ 1. Finally, substituting 
this in equation (^4|), we obtain 



d6n 6n 



dt pr'^ 



?/^^£^"P^-^°)^ 



.(c) 



Fk'> {e + i/2Y 

3uj 



(38) 



and thus the characteristic time-scale for the growth of angular speed is 

1 _ /• rf^ , .F'-:'(f. + l/2Y 



1 r d^ Ff{i+l/2) 
;^Ey^-^exp(-ro) 



(39) 



Figure 4 shows the growth time as a function of r. Note that the growth time near the 
top of the radiative zone is very short (of order a year). The growth time is inversely 
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proportional to the energy flux in waves, and therefore the time-scale for formation of the 
shear layer will be smaller if gravity waves are generated more efficiently than considered 
here (as should be the case if one considered also the waves excited by the convective 
overshooting); the magnitude of the differential rotation within the layer also increases with 
increase in the wave flux. 

The evolution of the angular speed at a given latitude below the convection zone is 
obtained by solving eq. ( p4|) and the results are shown in fig. 5. The initial rotation speed 
was taken to increase with depth below the convection zone, and the magnitude of the 
differential rotation over a layer of thickness O.OIR0 was considered to be ~ 0.05 nano-Hz 
initially (the mean rotation speed of the interior is 0.4 /iHz). We show results for several 
different co-latitudes (60" & 90°) and magnetic field strength (10^ Gauss & 8x10^ Gauss). 
The seed magnetic field we have considered lies in the horizontal plane along the ^-axis, 
which is the expected configuration for fields generated by differential rotation. The initial 
seed toriodal magnetic field we start with gets amplified with time as the magnitude of 
the shear in the layer just below the convection zone increases with time as a result of 
angular momentum deposit. However, to determine the time evolution of the magnetic field 
strength in a self-consistent way is outside the scope of this work. 

We see in fig. 5 that the time-scale for the formation of a strong shear layer is about 
20-30 years. The time-scale for the growth of the layer is a factor of about two smaller 
when the initial value of the differential rotation across the layer is 1 nano-Hz. Moreover, 
the growth time-scale is inversely proportional to the energy flux of waves. Considering the 
uncertainty in the initial 'seed' differential rotation rate and the energy flux in waves, the 
time-scale for the formation of the shear layer we flnd should be considered to be of the 
same order as the solar cycle period. 

In all of the cases shown in flg. 5, the rotation curve rises very sharply just below the 
convection zone, and this is followed by a more gently declining rotation curve. Below this 
layer, where positive angular momentum was deposited, lies another layer with negative 
differential rotation where retrograde gravity waves are absorbed. The sharp rise/decline of 
rotation in these two layers arises as a result of the rapid increase in wave damping rate 
with decreasing frequency in the local rest frame of the fluid i.e. an increase in the rotation 
rate causes prograde waves to be damped over a smaller distance scale which in turn 
increases the rotation rate and the feedback continues to increase the gradient of angular 
velocity until much of the angular momentum carried by prograde waves is absorbed. 
Since the frequencies of retrograde waves traveling inward in the star are increased as they 
propagate through a layer with negative gradient of differential rotation, they deposit their 
negative angular momentum deeper in the Sun leading to the formation of another thin 
and strong shear layer. The thickness of the shear layer is related to the damping length of 
low frequency gravity waves which carry bulk of the angular momentum flux. The distance 
from the top of the shear layer to the base of the convection zone is also about one damping 
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length. 

Shear layers at higher latitudes tend to be thinner and closer to the convection zone, 
and this is a result of the increase in the damping of wave^. Moreover, waves of large 
m and i do not propagate to high latitudes. Thus the formation of shear layer does not 
proceed at high latitudes. 

For waves of a fixed frequency, an increase in the strength of magnetic field generally 
leads to increase in the damping rate and, as a result, the shear layer is formed closer to the 
convection zone, is thinner, and grows more rapidly (see fig. 5). 

A self consistent calculation must consider the building up of the shear layer as well 
as the generation of magnetic field simultaneously. This is outside of the scope of this 
work. However, based on the physical understanding gained for the simpler problem 
considered here we expect that as the shear layer develops it twists poloidal magnetic field 
and generates toroidal magnetic field which in turn facilitates the damping of waves and 
prevents the redistribution of angular momentum in the layer by small scale turbulence 
(see §3.1). Both of these effects speed up the growth of angular momentum gradient in 
the layer. Another effect evident from figs. 5 & 6 is that the increase in the strength of 
the field causes the layer to move closer to the convection zone. When the magnetic field 
becomes sufficiently strong it prevents the propagation of low frequency waves, which are 
reflected at the top of the layer. These low frequency waves bounce back and forth between 
the convection zone and the shear layer and are absorbed after a few bounces thereby 
depositing their angular momentum and extending the shear layer upward closer to the 
bottom of the convection zone. A combination of shear instability and magnetic buoyancy 
could then lead to the ejection of the toroidal field from the shear layer and its deposition 
in the convection zone. 



3.1. The effect of shear instability on the buildup of shear layer 

So far, we have ignored the removal of angular momentum in the shear layer due to 
magnetic torque and the redistribution resulting from hydro dynamical instabilities. 

Magnetic torque acting on the shear layer removes angular momentum and therefore 
weakens the effectiveness of the gravity waves to build up the shear. The magnetic torque 
per unit volume Tb = r x (V x B) x B/47r ^ rBr{dB^/dr)/4TT. As long as the magnetic 
torque per unit area of the shear layer, rBrB^/Air, is small compared to the rate of angular 



''The spectrum of waves generated by turbulent convection has been taken to be independent of latitude 
which is not strictly valid and could modify this dependence. 
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momentum deposit in the layer we can ignore the effect of magnetic field. Since the angular 
momentum flux in prograde/retrograde waves is of order 2x10^^ g s~^, magnetic torque 
can be neglected so long as the product B^.B^ is less than about 5x10^ Gauss^, which we 
assume is the case. It should be noted that if gravity waves are excited by some process 
which is more efficient than the Reynold's stress considered in this paper, such as convective 
overshooting (cf. Fritts et al. 1998), then the limit on B^-B^ will increase correspondingly. 
The limit on BrB^ is, however, independent of the thickness of the shear layer since both 
the magnetic torque and the torque applied by gravity waves increase as the inverse of the 
thickness of the shear layer. 

The shear layer is expected to become unstable when the gradient of rotational velocity 
is sufficiently large, leading to the redistribution of angular momentum within the layer and 
possible mixing with the overlying convection zone. A shear layer in a stratified medium is 
unstable when the so-called Richardson criterion is satisfied (cf. Chandrasekhar 1961), 

g_dApl^ < i?,,, (40) 

p{dU/drf ^ ' 

where U is the horizontal velocity, Ap is the density difference between the perturbed fluid 
and the ambient medium, and Ric ~ | is the critical Richardson number for instability to 
set in. For an adiabatic perturbation, the above condition can be written as 

2 < R^c, (41) 



{dU/dz 

where A^^ = —g{dlnp/dr + g/c1) is the Brunt- Vaisala frequency (here, we will not consider 
the part of the Brunt- Vaisala frequency related to the mean molecular weight stratification 
as it is very small is the Sun's outer layers). Near the top of the radiative zone of the Sun 
N/2'K ^ 50/iHz. Thus a shear layer of thickness 5r becomes unstable when the differential 
rotation across this layer 5VL ^ 2N6r/rc, or when 6Q/2tt ^ 0.1/iHz for 6r/rc of 10^^. 

Even when the Richardson criterion for instability as described above is not satisfied, 
instability can still set in on small length scales due to the weakening of the buoyancy force 
as a result of thermal diffusion. The first elements to become turbulent are of size such that 
the time scale for thermal exchange with the surrounding fluid is small and the stabilizing 
effect of positive entropy gradient is reduced (see Zahn 1992, Maeder 1995). The modified 
Richardson criterion which includes heat diffusion is as follows 

/ 7] 



v + i 




N'<mc{^\ , (42) 



where r] = vlt/QK, v is the velocity of turbulent elements. It is their size, K is the thermal 
diffusivity defined in eq. (32), k, is opacity per unit mass, and Cp is specific heat per unit 
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inass, and a is the Stefan-Boltzmann constant. The turbulent diffusivity {D^) is given by 
those eddies which have the largest value of vlt and which satisfy (0), and it takes the 
following simple form: 

vlt „ .MU/drf 
D, ^ ^ ~ 2Ri,K^-^^ (43) 

(cf. Maeder 1995). This expression remains valid provided A^^/(df//dr) < ffi^. 

We include the effect of turbulent diffusivity on the evolution of the shear layer and 
the results are shown in fig. 6. Note that the main effect of the turbulence is to reduce the 
sharp velocity gradient at the top of the shear layer and to increase the thickness of the 
layer somewhat (compare figs. 5 & 6)Q. 

When the velocity gradient in the shear layer reaches the critical value so that 
\dU/dr\ = NRil^"^, the nature of the instability changes. The thermal diffusion is no longer 
required to destabilize the fluid and a dynamical instability sets in whose growth rate is 
of order of the rotation period (see e.g. Mac Donald 1983). The shear layer, which lies 
just below the convection zone, is likely to merge with the convection zone depositing the 
magnetic field that was build up in the layer as a result of differential rotation. This could be 
the start of a new magnetic cycle. That situation was never reached in our calculations since 
the diffusive shear instability was always large enough to redistribute angular momentum 
before condition (^) was satisfied. However, due to the lack of a reliable model, we have 
ignored the stabilizing effect of magnetic fields on the shear instability. We expect that 
magnetic fields will reduce the value for the turbulent diffusion coefficients we have used 
here, and therefore radial differential rotation will build up to larger values than shown in 
fig. 6 resulting in a more efficient generation of a toroidal magnetic field. In this case the 
dynamical instability condition is likely to be satisfied. The shear layer is also more likely 
to become dynamically unstable if the angular momentum fiux in gravity waves is larger 
than estimated in §2.3. 

The time-scale for the buildup of the shear layer and its distance from the base of 
the convection zone decreases with increasing latitude. Thus, dynamical instability and 
mixing of layer is expected to first occur at high latitudes and later in the equatorial region. 
However, a complete description of the latitude dependence would require a truly two 
dimensional model. Indeed, as there is no stabilizing stratification along horizontal layers, 
horizontal shear instabilities are much stronger than their vertical analogs, and we expect 
important horizontal redistribution of angular momentum to occur. 



'"^In this work, meridional circulation was ignored since it is known that, when differential rotation is 
strong, the momentum transport is dominated by shear turbulence (cf. Talon et al. 1997). 
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4. Discussion 

The gravity waves generated by turbulent convection are potentially very efficient in 
redistributing angular momentum in the radiative interior of the Sun. Most of the angular 
momentum flux is carried by the lowest frequency waves, which have periods of the order of 
the convective time-scale at the bottom of the convection zone (or ~ 0.3;uHz for the present 
day Sun). The total energy flux in these low period gravity waves is of the order of 0.03% of 
the solar luminosity and the angular momentum flux (Fl) in prograde waves is ~ 2 x 10^^ 
g s~^. We assume that prograde waves, which carry positive angular momentum flux, and 
retrograde gravity waves, which have negative angular momentum flux, are excited to the 
same amplitude and thus, the total net angular momentum flux at the top of the radiative 
interior is zero; in this case, gravity waves alone merely redistribute angular momentum in 
the interior. 

Kumar & Quataert (1997) and Zahn, Talon & Matias (1997) suggested that gravity 
waves generated by turbulent convection could bring the radiative interior of the Sun 
in corotation with the convection zone. The work presented here corrects a sign error, 
and adds one crucial element to their picture that was left out i.e. the formation of a 
double shear layer below the convection zone fllters out either prograde or retrograde waves 
depending on the sign of the initial velocity gradient below the convection zone. 

The thickness of these shear layers is of the order of the wave damping length (which 
is a few percent of the pressure scale height). The magnitude of the differential rotation 
(|(5i7|/27r) in the shear layers reaches about 0.4 /iHz when turbulent diffusion of angular 
momentum is not included. The magnitude of 6Q/2tt can be estimated using the relation 
|5n| ~ tFi/{pr'^5r) (this is applicable when waves with one sign of angular momentum are 
absorbed in the layer), and we flnd it to be 0.1 /xHz after a period of a few years when 
the thickness of the shear layer is 3 x lO~^i?0. The thickness is uncertain by a factor of 
a few due to the uncertainty in spectrum of gravity waves and the angular momentum 
redistribution as a result of shear and magnetic instability. 

As the angular velocity increases in the layer, it becomes unstable to small length 
scale perturbations which cause mixing of the layer, and thus contributes to the mixing of 
elements near the top of the radiative interior of the Sun as seems to be required by the 
helioseismic inversion of the sound speed (Basu, 1997). 

If only "diffusive" turbulence is present (as in our calculations), the flrst shear layer 
lying at the top of the radiative zone will eventually merge with the surface convection 
zone and the second shear layer, of opposite sign, will replace it at the top of the radiative 
zone, reversing the structure, much as in the biennial oscillation, which is well known in 
atmospheric sciences (cf. Holton & Lindzen, 1972). Taking into account the stabilizing 
effect of a toroidal magnetic fleld could change this situation as stronger shears would 
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result, leading perhaps to the appearance of a more drastic "dynamical" instability which 
might destroy much of the double shear layer. 

Most of the angular momentum deposited in these shear layers is transferred back 
into the convection zone once they become turbulent, but further work is needed to check 
whether there will be a net amount of angular momentum taken out of the radiation zone, 
and how it will depend on the angular velocity gradient. This issue must also be settled in 
order to get a clear picture of the interaction between these short-lived shear layers and 
the quasi-stationary tachocline, which is the somewhat thicker transition layer that links 
the differential rotation of the convection zone to the almost uniform rotation below (see 
Spiegel & Zahn 1992). 

A poloidal magnetic field threading the shear layer will get twisted as the differential 
rotation in the layer increases and, as a result, toroidal magnetic field will be generated. 
The increasing toroidal field in turn accelerates the rate of angular momentum deposit in 
the layer, and inhibits small length scale shear instability which reduces the rate of growth 
of the angular momentum gradient. Once the magnetic field in the shear layer becomes 
strong enough to reflect low frequency or high azimuthal wavenumber waves, these waves 
travel back and forth between the shear layer and the convection zone and their angular 
momentum is absorbed in this intermediate layer in a few transit time. The result of this 
is to extend the shear layer closer to the convection zone. At some stage in its evolution, 
the shear layer becomes unstable and mixes with the convection zone and the magnetic 
field is the layer floats upward to start a new magnetic (half) cycle. Once the layer has 
relaxed and settled back to a stable configuration, after having given up much of the 
angular momentum and magnetic field that it had contained, the whole process starts all 
over again. Let us point out that this model predicts cyclic reversals of the toroidal field 
with a steady poloidal field at the base of the convection zone, which could well be a fossil 
field, whereas in classical dynamo theory the poloidal field too undergoes such reversals. 
We note also that waves of large azimuthal order (m) do not reach high latitude. Therefore, 
the formation of the shear layer and the generation of magnetic field in not likely to occur 
at high latitudes which is consistent with the lack of sunspots there. 

Incorporating these various processes in a self consistent manner in the evolution 
calculation is, however, very complicated and has not been considered here. But it is 
tempting to link the cyclic build-up of these shear layers, which produce toroidal field of 
alternating sign, with the solar magnetic cycle, since the time-scales are of the same order. 
Future helioseismic observations and inversions with improved resolution should be able to 
determine the presence of a double shear layer at the top of the radiative interior predicted 
here. 
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Figure Captions 

FIG. 1. — The radial wavenumber {krRo) as a function of wave frequency for B = 10^, & 
8x10^ Gauss, and for latitude {le)= 0° k 30". The sohd line is for {B = 10^, Ig = 0°), the 
dotted hue {B = 10^ k = 30"), dashed line for {B = 8 x 10^ Ig = 0°), and the dash-dot hue 
is for {B = 8 X 10^, Ig = 30°). In all of these cases the horizontal wave number was taken to 
be kh = 50/Rq, m = 30, Q/2n = 0.4/iHz, p = 0.2 g cm~^, and the Brunt- Vaisala frequency 
N = 100/iHz (the last three values correspond to the top of the solar radiative zone). 

FIG 2. — The left panel shows energy luminosity per unit radian frequency in gravity waves 
of degree £ = 1 as a function of frequency (energy luminosity is defined to be the product 
of energy flux and the surface area). The right panel shows the frequency integrated wave 
action, or angular momentum flux divided by m, as a function of i; the frequency integration 
was carried out over the interval 0.6-3.0 /iHz (most of the contribution to the wave action 
comes from frequency close to 0.6 /iHz, see the left panel). We have used the solar model 
of J. Christensen-Dalsgaard in this as well as all the other numerical calculations shown in 
this paper. 

FIG 3. — Damping length for gravity waves divided by the wavelength of the wave is shown 
as a function of wave frequency for i = 2 (solid line), 15 (dotted line), and 50 (dashed line). 

FIG 4. — The time-scale for exponential growth of angular velocity gradient below the 
solar convection zone as a function of r for two different values of the magnetic field 
B = 10^, & 8x10'^ Gauss, and latitude {Ig =) 0° & 30°. The top of the radiative zone is 
at t/Rq = 0.71195. The growth time-scale was calculated by integrating over frequency 
range of 0.6-3.0 /iHz and the maximum i was taken to be 55. The solid line is for 
{B = 10^, Ig = 0°), the dotted line {B = 10^, Ig = 30°), dashed hue for {B = 8 x 10^, Ig = 0°), 
and the dash-dot line is for (B = 8 x 10^ Ig = 30°). 

FIG 5. — Evolution of the shear layer at two different latitudes (0° and 30°) and for B = 10^ 
and 8x10^ Gauss. The initial rotation profile is taken to increase with depth and the 
magnitude of the initial differential rotation is about 2 nano-Hz over the range of r shown 
in these plots. The thin solid line is the rotation profile after 10 yr, dotted line at 20 yr, 
dashed line at 30 yr and dash-dot line at 40yr. Turbulent diffusion of angular momentum 
in the shear layer is not included in these calculations. 

FIG 6. — Evolution of the shear layer. The top left pannel shows the difference between 
a calculation including turbulent diffusion (thin continuous line) and excluding turbulent 



This preprint was prepared with the AAS h^TJrjX macros v4.0. 



-21- 



diffusion (thick dashed hne) at the equator and for B = 10^ Gauss at 20, 30 and 40 years. 
The other three pannel show the shear layer including turbulent diffusion at 0° and 30° and 
B = 10^ and 8x10^ Gauss. The initial rotation profile is as in Fig. 5. Ages are indicated on 
the plot. 
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